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ABSTRACT 

We discuss Spherical Needlets and their properties. Needlets are a form of spherical 
wavelets which do not rely on any kind of tangent plane approximation and enjoy good 
localization properties in both pixel and harmonic space; moreover needlets coefficients 
are asymptotically uncorrelated at any fixed angular distance, which makes their use in 
statistical procedures very promising. In view of these properties, we believe needlets 
may turn out to be especially useful in the analysis of Cosmic Microwave Background 
(CMB) data on the incomplete sky, as well as of other cosmological observations. As 
a final advantage, we stress that the implementation of needlets is computationally 
very convenient and may rely completely on standard data analysis packages such as 
HEALPix. 
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1 INTRODUCTION 

Over the last few years, wavelets have emerged as one 
of the most powerful tools of CMB data analysis, find- 
ing applications in virtually all areas where statistical 
methods are required; a very incomplete list of references 
should i nclude testing fo r non- Gaussianity (see lVielva et al.l 
1 2004): ICabella et al] l|20040 . foreground subtraction 
( Hansen et al. 1 2006t )). point sour ce detection JSanz et al.l 
( 2006t)). componen t separation (|Moudden et aL J2005I); 
Starc k et al.l ((2006)), polarization analysis I Cabella et al.l 
|2007l )) and many others. The reason for such a strong in- 
terest is easily understood. As it is well-known, CMB mod- 
els are best analyzed in the frequency domain, where the 
behaviour at different multipoles can be investigated sepa- 
rately; on the other hand, partial sky coverage and other 
missing observations make the evaluation of exact spherical 
harmonic transforms troublesome. The combination of these 
two features makes the time-frequency localization proper- 
ties of wavelets most valuable. 

Despite the wide agreement on their importance as 
a data analysis instrument, the derivation of an optimal 
wavelets basis on the sphere is still an open issue for re- 
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search. Many efforts have been undertaken in this area, 
most of th em being based upon the so-calle d tangent plane 
approach (jAntoine &: Vanderghevnstl fl9 99V In this frame- 
work, a flat sky approximation is entertained locally, and 
then some form of standard plane wavelets ar e implemented. 
Directional wavelets have been advocated bv lMcEwen et all 
(2006, 2007), again by means of a tangent plane approxima- 
tion. An interesting atte mpt to overcome the tangent plane 
approximation is due to ISanz et all l|2006l ). 

A new approach to spheri cal wavelets was in troduced 
in the statistical literature by iBaldi et al.l l|2006l ). adapt- 
ing tools propos e d in the functional analysis literature by 
iNarcowich et all (|2 006k th e first application t o CMB data 



is due to 



Pietrobon et al.l (|2006l ) , (see also iGuilloux et al.l 



l|2007h : lBaldi et all jjooj )). The idea is to focus on so-called 
needlets, to be described in the following section. Needlets 
enjoy several features which are not in general granted by 
other spherical wavelets construction; we anticipate some of 
these features, which we shall investigate more deeply in the 
Sections to come. More precisely: 

a) they d o not rely on an y tangent plane approxima- 
tion (compare ISanz et al.ll2006l ). and take advantage of the 
manifold structure of the sphere; 

b) being defined in harmonic space, they are computa- 
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tionally very convenient, and natively adapted t o standard 
packages such as HEALPbQ (|G6rski et alj|2005h : 

c) they allow for a simple reconstruction formula 
(see Eq. [5}, where the same needlets functions appear 
both in the direct and the inverse transform (see also 
IKerkvacharian et al, I l2007l ). This property is the same as for 
spherical harmonics but it is not shared by other wavelets 
systems such as the well-known Spherical Mexican Hat 
Wavelets (hereafter SMHW); 

d) they are quasi-exponentially (i.e. faster than any 
polynomial) concentrated in pixel space, see Eq. [6] below; 

e) they are exactly localized on a finite number of mul- 
tipoles; the width of this support is explicitly known and 
can be specified as an input parameter (see Eq. [TJ ; 

f) random needlets coefficients can be shown to be 
asymptotically uncorrelated (and hence, in the Gaussian 
case, independent) at any fixed angular distance, when the 
frequency increases. This capital property can be exploited 
in several statistical procedures, as it allows to treat needlets 
coefficients as a sample of independent and identically dis- 
tributed coefficients on small scales, at least under the Gaus- 
sianity assumption. 

The aim of this paper is to discuss more thoroughly the 
implementation of needlets, compare it with other wavelets 
(namely, the SMHW) and investigate their properties by 
means of Monte Carlo simulations. In Section [2] we describe 
the numerical implementation, taking care to discuss the fea- 
tures of the construction that ensure the above-mentioned 
properties. In Section [3] we discuss the relationship between 
the localization properties in frequency and pixel spaces; we 
also discuss the trade-off between the two, which from the 
operational point of view relates to the issue of an opti- 
mal choice of a so-called "bandwidth" parameter B. Some 
comparisons are made with existing techniques to deal with 
partial sky coverage, and more precisely with tophat bin- 
ning procedures and SMHW. Some discussion on possible 
applications and directions for further research is collected 
in Section [4] 



2 THE NUMERICAL IMPLEMENTATION OF 
NEEDLETS 

We start by outlining briefly the c onstruction of a needlet s 
basis. More details can b e found in iNarcowich et all <|2006h . 
and in lBaldi et al. (2006). We shall discuss the details of the 
construction step by step, in order to provide to potential 
users a clear recipe for needlets implementation. 

We first recall that the spherical needlet (function) is 
defined as 



Needlet shape: B = 3, j = 8 



(1) 



here, we use {£jk} to denote a set of cubature points on the 
sphere, corresponding to frequency j; in practice, we shall 
identify these points with the pixel centres in HEALPix. 
Also, Xjk denotes the cubature weights, which for simplicity 
can be envisaged as 1/N V , N p denoting the number of pixels 
fsee lPietrobon et aij|2009 ). 
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Figure 1. Needlets in pixel space. B = 2, j 



Intuitively, needlets should be viewed as a convolution 
of the projection operator J2m=-e Ytm(,7)Yi m (£jk) with a 
suitably chosen window function &(•) Special properties of 
&(■) ensure that the needlets enjoy quasi-exponential local- 
izatio n properties in pixel space. Formally, w e must ensure 
that l|Narcowich et al.ll2006l ; iBaldi et alj|2006h : 

(i) 6 2 (-) has support in [^,B], and hence b{-~) has sup- 
port in £ <E [B^^B^ 1 ] 

(ii) the function &(•) is infinitely differentiable in (0, 00). 

(iii) we have 



1 for all £ > B. 



(2) 



It is immediate to see that property (i) ensures the 
needlets have bounded support in the harmonic domain; 
property (ii) is the crucial element in the derivation of the 
localization properties, which we shall illustrate in the fol- 
lowing section. Finally, property (iii) is necessary to estab- 
lish the reconstruction formula which we shall discuss below; 
functions such as b 2 (-) are called partitions of unity. 

There are of course many possible constructions satis- 
fying (i-iii); indeed an interesting theme for future research 
is the derivation of optimal windows sa tisfying these three 
conditions fcompare lGuilloux et al.ll2007l ). An explicit recipe 
for the construction of &(•) is as follows. 

STEP 1: Construct the function 



/(*) = 



exp(- 
0, 



-1 < t < 1 
otherwise 



It is immediate to check that the function /(•) is C°° and 
compactly supported in the interval (—1,1) 
STEP 2: Construct the function 



JV(«)<ft 
£ /(*)*' 



The function ?/>(■) is again C°°; it is moreover non-decreasing 
and normalized so that 1) = , = 1 
STEP 3: Construct the function 
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<p(t) = { va- 
o 



_2B 

D 



if 
if 
if 



«S 



1_ 

B 
«S 1 



t > 1 



Here we are simply implementing a change of variable so 
that the resulting function ip(-) is constant on (0, B~ ) and 
monotonically decreasing to zero in the interval (_B _1 ,1). 
Indeed it can be checked that 



identity in Eq.|4]has indeed been verified by means of numer- 
ical simulations and implicitly provides the correct normal- 
ization for needlets. It is basically a consequence of the pe- 
culiar partition-of-unity property of &(•) of Eq. [2] Of course 
this property is not generally shared by other constructions 
such as SMHW, where the wavelets functions are normal- 
ized to unity in the real domain. Eq. [4] is related to a much 
more fundamental result, i.e. the reconstruction formula 



2B 



B 



-it- 



1= 



1 

-1 



for 
for 



t = 1 



and 



m B 



<p{i) = ^(-i) = o 

STEP 4: Construct 
b 2 (0 = ¥>(§) - <p{Q 

The expression for b 2 (-) is meant to ensure that the function 
satisfies the partition-of-unity property of Eq. [21 Needless to 
say, for b(£) = {<p(-g) — ¥>(£)} ^ we take the positive root. 



Random needlets coefficients are hence given by 



A ^E 6 ^ E {J g2 T(VYi m tf)df^Ye m (i 



(3) 



It is very important to stress that, although the needlets 
do not make up an orthonormal basis for square integrable 
functions on the sphere, they do represent a tight frame. In 
general, a tight frame on the sphere is a countable set of 
functions {e^} such that, for all square integrable functions 
on the sphere / £ L 2 (S 2 ), we have 



E 



so that the norm is preserved. Of course, this norm- 
preserving property is shared by all orthonormal systems; 
however, frames do not in general make up a basis, as 
they admit redundant elements. They can be viewed as 
the closer system t o a ba si s, for a given redu ndancy, see 
iHernandez fc Weiss! l|l99rj ). iBaldi et all (|2006T ) for further 
definitions and discussion. In our framework, the norm- 
preserving property becomes 



j,k f=l 

where 



2£ + 1 

47T 



(4) 



Ct 



21 + 



—y 

+i ^ 



(I < n 



(•») 



which in turn is a non-trivial consequence of the careful con- 
struction leading to Eq. [2] As mentioned before, the simple 
reconstruction formula of Eq. is typical of tight frames 
but does not hold in general for other wavelets systems. It 
is easy to envisage many possible applications of this result 
when handling masked data. 



3 PROPERTIES AND COMPARISONS 

The following qua si-exponential lo c alizat ion property of 
needlets is due to iNarcowich et al.1 (|2006T ) and motivates 
their name. 

For any k = 1,2, ... there exists a positive constant Ck 
such that for any point 7 e S 2 we have 



l^fc(7) I 



Ck 



B J 



(1 + Biarccos((7,00) t 



(6) 



We recall that arccos((7, f jk)) is just the natural dis- 
tance on the unit sphere between the points (7, The 
meaning of Eq. [6] is then clear: for any fixed angular dis- 
tance, the value of ijijkif)) goes to zero quasi-exponentially 
in the parameter B. The resulting trade-off in the behaviour 
over the harmonic and real spaces is expected: smaller values 
of B correspond to a tighter localization in harmonic space 
(less multipoles entering into any needlet), whereas larger 
values ensure a faster decay in real space. 

Due to their localization properties, needlets are espe- 
cially useful in the analysis of partial sky coverage. In fact, 
in view of Eq. [S]we expect the value of needlets coefficients 
to be mildly affected by the presence of gaps in the maps. 
To illustrate this important feature, we plot the quantity 



Xjk 



jk,mask 



- M 2 ) 



{flu) 



(7) 



in Fig. [5J where the KpO mas that is used to remove 
Galactic emission and point sources from WMAP data 
(roughly 75% of the sky), is applied. The expected values 
of Eq. [7] are again evaluated by means of 100 Monte Carlo 
simulations; in particular we focus on needlets coefficients 
corresponding to B — 2.72 and j — 5, which amounts to a 
range in multipoles space in the order of 58 < £ < 398. To 
put our results in perspective, in the same Figure we show 
analogous findings with the use of a tophat binning filter 
and SMHW. We remind the SMHW formula 



1 



2nN(R) 



1 



is the raw angular power spectrum of the map Tiff). The 
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4 D. Marinucci et al. 




KpO mask effect on tophat binning coefficients 




KpO mask effect, on SMHW coefficients - R - 28 aremin 




fitL^ » **** jS^JZ, f *» 




Figure 2. Effect of incomplete sky coverage, modeled by the 
WMAP KpO mask, visualized by plotting on a sky map the quan- 
tity defined in Eq. [7] From top to bottom, the result for needlets, 
flat binning, and SMHW (28')- 



- (!) 2 eX p(-2/ 2 /2-R 2 ; 

where y — 2 tan (9/2 (8 is the polar angle), R is the 
scale of convolution and N( R) a normalization factor 
|Martmez-Gonzalez et al.ll2002n . 

Under these circumstances, needlets coefficients are well 
localized, but slightly sensitive to the mask. Indeed, only 
56% of the pixel are changed by less than 0.1; SMHW co- 
efficients perform a bit better (73%) while a simple tophat 
binning fails completely (only 6%). The difference between 
the two wavelets bases can be due to the different power that 
they give to multipoles (see Fig. [4]). In fact the performance 
of needlets can be improved choosing the appropriate B J , 
that defines the optimal shape for the window &(•), given the 



multip oles range of interest. For details see iGuilloux et all 
((2003). In the same paper, the authors argue that an op- 
timal filter can be adapted to deconvolve a specific mask: 
this property provides a further degree of flexibility to the 
needlet s toolbox. 

In iBaldi et al] (|2006l ). another relevant property of 
needlets coefficients was discussed, namely their asymptotic 
uncorrelation at any fixed angular distance, for growing fre- 
quencies j. More explicitly, at high frequency needlets coef- 
ficients can be approximated as a sample of identically dis- 
tributed and independent (under Gaussianity) coefficients, 
and this property opens the way to a huge toolbox of sta- 
tistical procedures for CMB data analysis (for instance, for 
testing Gaussianity and isotropy). Also, in view of Eq. [3] for 
full sky maps and in the absence of any mask we should ex- 
pect the theoretical correlation to be identically zero when- 
ever |ji -j 2 \ > 2. 

Let us define the realized correlation between two dif- 
ferent scales ji, j2 as 



Piii 



(8) 



By using Eq. [3] one has that: 

(Phki3 j2k ) = b (£/B jl ) b{£/B j2 ) G t K "*> 



where Ci is the underlying CMB angular power spectrum 
and the coupling kernel Ktfr nm.' is defined in t erms of the 
observed mask W(7) (see e.g. iHivon et al] |2002): 



Note that, in the absence of gaps (i.e. W(j) = 1), 
y"] m Kttmm reduces simply to {21 + 1). 

Equation (JHJ can not be expected to be reproduced ex- 
actly, due to numerical approximations; in particular, we 
should stress that theoretical results are derived under the 
assumption that needlets coefficients are evaluated at ex- 
act cubature points, so that the {a; m } are precisely recon- 
structed from the maps. Of course, this is not the case in 
practice; however, we do expect small and vanishing values 
for ji <C J2- At the same time, we expect this correlation 
to increase on the presence of sky cuts, but less so than for 
other bases. Here, we want to illustrate the practical rele- 
vance of this mathematical results by means of simulations 
on the correlation coefficient. More precisely, we computed 
the quantity © by performing a Monte Carlo over 100 sim- 
ulations. Our findings are shown in Tables [TJ [2] 

We view these results as very encouraging. In the ab- 
sence of a mask, the correlation coefficient is by any practical 
means virtually negligible for all frequency distances greater 
or equal than 2, while at distance Aj = 1 the correlation is 
around ~ 0.25 in good agreement with eqn. [5] which pre- 
dicts 0.22 for our input parameters. In the presence of sky 
cuts, the performance deteriorates as expected only at low 
j where it exceeds a few percentage points, as shown for our 
simulations in the case of the KpO mask. A computation 
analogous to ((8]) yields for SMHW the theoretical results 
reported in Table [3] note how we have non zero values at 
all lags. Numerical results to support the theoretical find- 
ings are provided by Tables |4j [5] We believe these compared 
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Table 1. Needlets correlation parameter. B = 2.72 without gaps 



Hi' 1 2 3 4 5 

1 1.000 0.275 0.001 0.001 0.003 

2 - 1.000 0.248 0.001 0.001 

3 - - 1.000 0.268 0.001 

4 - - - 1.000 0.242 

5 - - - - 1.000 



Table 2. Needlets correlation parameter. B = 2.72 with gaps 

Hi' 1 2 3 4 5 

1 1.000 0.420 0.140 0.040 0.060 

2 - 1.000 0.335 0.023 0.001 

3 - - 1.000 0.291 0.004 

4 - - - 1.000 0.252 

5 - - - - 1.000 



results strongly support the potential of needlets for the im- 
plementations of statistical procedures, where uncorrelation 
properties are clearly a very valuable asset. 

As a further comparison, we evaluated the domains in 
harmonic and real spaces for needlets, tophat binning and 
SMHW. In particular we normalized the three bases to have 
roughly an equal area in the harmonic domain, paying at- 
tention to have the maximum of the power in a similar range 
of multipoles. Results are plotted in Fig. U [S] It is evident 
how SMHW and needlets outperform tophat binning by two 
order of magnitudes in terms of localization in real space: in- 
deed in this domain the two wavelets constructions perform 
quite similarly. Moreover, in Figure [7] we computed the an- 
gle where the integral of the filter functions in pixel space 
reaches 68%, 95% and 99% of the total area, respectively. 
Again, it is immediate to check how at every scale needlets 
outperform very clearly a simple binning approach; on the 
other hand, SMHW seems slightly more concentrated in this 
setting. The linear trend for needlets in the log-log plot is 
a direct conseguence of their construction, and in particular 
of the functional dependence on £/B J . 

On the other hand the advantage of needlets over 
SMHW emerges quite clearly in the harmonic domain. More 
precisely, after normalizing the two methods to be centred at 
the same angular scale, with roughly the same total power, 
the needlets support seems clearly more concentrated than 
SMHW. In particular we stress how SMHW suffer from 

Table 3. Theoretical correlation (full sky) for needlets and 
SMHW 



corr / Aj 





1 


2 


3 


i 


Needlets 


1.000 


0.220 


0.000 


0.000 


0.000 


SMHW 


1.000 


0.500 


0.100 


0.014 


0.002 



Window functions in harmonic domain 
1.2 i i 




500 



Figure 4. The red solid line represents needlets window function, 
b{-~) in harmonic space for B = 2.72, j = 5. The blue dashed 
and green dot-dashed lines provide the tophat and the SMHW 
window functions, respectively. The SMHW corresponds to a scale 
R = 28/ in pixel space. 

Filter profiles in real space 




-2-10 1 2 



a 

Figure 5. Behaviour of needlets (solid red), SMHW (dot-dashed 
green) and tophat binning (blue dashed) in pixel space. The angle 
in horizontal axis is measured in radiants. 

"leakage" by the very low multipoles, i.e. exactly those most 
affected by sky cuts and cosmic variance. No such leakage 
occurs for needlets. 



4 CONCLUDING REMARKS 

In this paper, we have described a new construction for a 
spherical wavelets frame, and we have illustrated several of 
its properties; most notably localization properties in the 
real and harmonic spaces, uncorrelations of the resulting 
random coefficients, and independence from any tangent 
plane approximations. Moreover, needlets enjoy a direct re- 
construction property which allows analysis and synthesis 
to be implemented directly and in a computationally conve- 
nient manner. Each of these properties has been illustrated 
by means of Monte Carlo simulations. The encouraging re- 
sults reported suggest that needlets can become a valuable 
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Table 4. SMHW correlation parameter without gaps. The scale R is given in arcmin. 



R/R' 


1792 


896 


448 


224 


112 


56 


28 


1792 


1.000 


0.503 


0.109 


0.016 


0.002 


0.0002 


0.00003 


896 




1.000 


0.500 


0.099 


0.014 


0.002 


0.0002 


448 






1.000 


0.510 


0.103 


0.014 


0.002 


224 








1.000 


0.511 


0.104 


0.014 


112 










1.000 


0.513 


0.107 


56 












1.000 


0.519 


28 














1.000 



Table 5. SMHW correlation parameter with gaps. The scale R is given in arcmin. 



R/R' 


1792 


896 


148 


224 


112 


56 


28 


1792 


1.000 


0.496 


0.113 


0.022 


0.005 


0.002 


0.0007 


896 




1.000 


0.520 


0.115 


0.021 


0.005 


0.002 


448 






1.000 


0.523 


0.114 


0.021 


0.005 


224 








1.000 


0.520 


0.114 


0.020 


112 










1.000 


0.522 


0.116 


56 












1.000 


0.527 


28 














1.000 



Filter profiles in real space - details 



Localization comparisons 
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CD 

s 



10.0 r 




100 



Figure 6. Here in this figure we provide details of the behaviour 
in pixel space over the relevant range, i.e. the region where the 
three functions exceed 0.001. Lines have the same meaning as in 
the previous figures. 



tool in the several areas of CMB data analysis where other 
wavelets have already proved useful. 
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Figure 7. The three lines represent the angle at which the area 
of needlet, in red, or the tophat, in blue, filter reaches the 99%, 
95%, 68% of the total area as a function of the peak multipolc 
in each window function. The latter corresponds directely to a 
given j for needlets and to the scale R for SMHW; for the tophat 
window the central i in the band is taken. 
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